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1. Introduction 

In a recent issue of this Monthly, Hetzel, Liew, and Morrison [4] pose a rather natural ques- 
tion: what is the probability that a random n x n integer matrix is diagonalizable over the field of 

rational numbers? 

Since there is no uniform probability distribution on Z, we need to exercise some care in inter- 
preting this question. Specifically, for an integer k >1, let Jfc = {—k, — /c + 1, . . . , A; — 1, A;} be 
the set of integers with absolute value at most k. Since is finite, we are free to choose each entry 
of an n X n matrix M independently and uniformly at random from X^, with each value having 
probability l/{2k + 1). The probability that M has a given property, such as being diagonalizable 
over Q, is then a function of the parameter k; we consider how this function behaves as A; — > oo. 
In particular, if the probability converges to some limiting value then it is natural to think of this 
limit as the "probability" that a random integer matrix has that property. 

We refer the reader to the article of Hetzel et al for an interesting discussion of some of the 
issues raised by this interpretation of probability as a limit of finitary probabilities (in doing so we 
lose countable additivity and hence the measure-theoretic foundation of modem probability theory 
after Kolmogorov). From a pragmatic viewpoint, this cost is outweighed by the fact that many 
beautiful number-theoretic results are most naturally phrased in the language of probability: for 
instance, the celebrated Erdos-Kac theorem [1] states that the number of prime factors of a positive 
integer n behaves (in an appropriate limiting sense) as a normally-distributed random variable with 
mean and variance both equal to log log n. (In this article we always mean the natural logarithm 
when we write log.) 

For any given integers n>2 and k >1, the set of random nxn matrices with entries in X^ is a 
finite probability space; it will be convenient to compute probabilities simply by counting matrices, 
so we introduce some notation for them. Let M.n{k) denote the set of all n x n matrices whose 
entries are all in X^; then we are choosing matrices uniformly from M.n{k), which has cardinality 
exactly {2k + 1)" . The probability that a random matrix in M.n{k) satisfies a particular property 
is simply the number of matrices in M.n{k) with that property divided by {2k + 1)" . 

For a given integer A, let A4^(A;) denote the set of all matrices in Mn{k) that have A as an 
eigenvalue. Note that in particular, (A;) is the subset of singular matrices in M.n{k). Likewise, 
we denote set of the matrices in A4.n{k) having at least one integer eigenvalue hy A4^{k) = 
Uagz -^n(^)- The probability that a random matrix in M.n{k) has an integer eigenvalue is thus 
\Ml{k)\/{2k + l)-\ 

Our main result affirms and strengthens a conjecture made in [4]: /or any n > 2, the probability 
that a random n x n integer matrix has even a single integer eigenvalue is 0. We furthermore give 
a quantitative upper bound on the decay rate of the probability as k increases. It will be extremely 
convenient to use "Vinogradov's notation" to express this decay rate: we write f{k) <C g{k) if 
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there exists a constant C > such that |/(A;)| < Cg{k) for all values of k under consideration. 
Notice for example that if fiik) < g{k) and /2(A;) < g{k), then f\{k) + /2(A';) < ^((A;) as well. 
If this constant can depend upon some auxiliary parameter such as e, then we write f{k) -C^ g{k); 
for example, for A; > 1 it is true that log k k^ for every e > 0. 

Theorem. Given any integer n > 2 and any real number e > 0, the probability that a randomly 
chosen matrix in A4n{k) has an integer eigenvalue is <tin,s ^/k^~^- In particular, the probability 
that a randomly chosen matrix in Aln(^) diagonalizable over the rational numbers is <^n,e 

Given an integer matrix M e Mn{k),?L necessary condition for it to be diagonalizable over Q is 
that all of its eigenvalues are rational. Moreover, since the characteristic polynomial det(A/ — M) 
is monic with integer coefficients, the familiar "rational roots theorem" [11, §4.3] implies that 
every rational eigenvalue of M must be an integer. Hence any matrix that is diagonalizable over 
the rationals must certainly belong to M^{k), and so the second assertion of the theorem follows 
immediately from the first. 

The special case n = 2 of the theorem was obtained in [4] and also earlier by Kowalsky [7]. 
Unravelling the ^-notation, the theorem states that there exists a constant C, possibly depending 
on n and e, such that \Ml{k)\/\Mn{k)\ < C/k^-' for all A; > 1. Note that \Mn{k:)\ A;"' 
(with the implied constant highly dependent on n), and so the theorem also gives an upper bound 
for the number of matrices 'mM.n{k) with at least one integer eigenvalue, namely 

The key tool used to establish the theorem is the following related estimate for the number of 
singular matrices in Mn{k): 

Lemma 1. Given any integer n > 2 and any real number £ > 0, the probability that a random 
matrix in Ain{k) is singular is <^n,e 1/A;^~^. In other words, \A4^{k) \ <^n,e ^" 

This upper bound is essentially best possible when n — 2: we can show that the cardinality of 
Ai2{k) is asymptotic to (QG/tt^) A;^ log k, so that the probability that a matrix in M.2{k) is singular 
is asymptotic to (6 log k)/n'^k'^. (We discuss the 2x2 case in more detail in Section 3.) However, 
we expect that the upper bound of Lemma 1 is not sharp for any n > 2. The probability that 
a matrix has a row consisting of all zeros, or that two of its rows are identical, is 1/A;" (up to 
a constant depending only on n), so the upper bound cannot be any smaller than this. It seems 
reasonable to conjecture that the true proportion of singular matrices in M.n{k) decays as 1/A;"~^, 
but for higher dimensions we do not know what the correct analogy should be of the precise rate 
(log k) I k"^ that holds in dimension 2. 

We remark that in this paper, we are considering the behaviour of \Ain{k)\/\Mn{k) \ when n is 
fixed and k increases. It is perfectly natural to ask how this same probability behaves when k is 
fixed and n increases, and there is a body of deep work on this problem. Komlds [6] proved in 1968 
that this probability converges to as n oo, answering a question of Erdos. In fact Komlos's 
result holds, not just for the uniform distribution on X^, but for an arbitrary non-degenerate distri- 
bution on R. (A degenerate distribution is constant almost surely, so the non-degeneracy condition 
is clearly necessary.) Slinko [14] later established a quantitative decay rate of <tik \fn. An ex- 
ponential decay rate of (0.999)" for the case of {±l}-matrices was estabUshed by Kahn, Komlos, 
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and Szemeredi [5] and improved to (3/4)" by Tao and Vu [15]. Very recently Rudelson and Ver- 
shynin [12] have established exponential decay for a very wide class of distributions, including the 
uniform distribution on X^. 

2. Determinants, singular matrices, and integer eigenvalues 

We begin by proving a lemma that we will use repeatedly in the proof of Lemma 1. It shows that 
the probability that a 2 x 2 matrix is singular remains small (that is, the probability is ^ A:^^+^ just 
as in Lemma 1, albeit with a different implied constant) even if we choose the entries randomly 
from arbitrary arithmetic progressions of the same length as Xk. 

Lemma 2. Fix positive real numbers a and e. Let k be a positive integer, and let Li{x), ^2(2;), 
L^lx), and I/4(x) be non-constant linear polynomials whose coefficients are integers that are <^a 
/c" in absolute value. Then the number of solutions to the equation 

L,{a)L2(b) ^ Ls(c)L4d), (1) 

with all of a, b, c, and d in X^, is <^a,e k^'^^. 

Proof First we consider the solutions for which both sides of equation (1) equal 0. In this case, at 
least two of the linear factors -Li(a), -^2(6), -^3(0), and L4{d) equal 0. If, for example, Li{a) = 
and i^3(c) = (the other cases are exactly the same), this completely determines the values of a 
and c; since there are 2/c + 1 choices for each of b and d, the total number of solutions for which 
both sides of equation (1) equal is <S k'^. 

Otherwise, fix any values for c and d for which the right-hand side of equation (1) is nonzero, a 
total of at most {2k + 1)^ <^ k"^ choices. Then the right-hand side is some nonzero integer that is 
at most (/c ■ A;" + A;")^ < 4^^"*"^" in absolute value, and Li(a) must be a divisor of that integer. 

It is a well-known lemma in analytic number theory (see for instance [9, p. 56]) that for any 
S > 0, the number of divisors of a nonzero integer n is <^s In particular, choosing S — 

e/{2 + 2a), the right-hand side of equation (1) has <^a,e ^/^f^2+2ay/{2+2a) divisors to serve 

as candidates for Li{a); each of these completely determines a possibility for a (which might not 
even be an integer). Then the possible values for .^2(6) and hence b are determined as well. We 
conclude that there are a total of <^a,e k"^^^ solutions to equation (1) as claimed. □ 

Remark. It is not important that the Lj be linear polynomials: the above proof works essentially 
without change for any four non-constant polynomials of bounded degree. We will not need such 
a generalization, however, as the determinant of a matrix depends only linearly on each matrix 
element. 

The next ingredient is a curious determinantal identity, which was classically known but at 
present appears to have fallen out of common knowledge. Before we can state this identity, we 
need to define some preliminary notation. For the remainder of this section, capital letters will 
denote matrices, boldface lowercase letters will denote column vectors, and regular lowercase 
letters will denote scalars. 

Let /„ denote the n x n identity matrix, and let denote the jth standard basis vector (that is, 
the jth column of /„). Let M be an n x n matrix, and let uij denote its jth column and rriij its ijth 
entry. Note that Mej = uij by the definition of matrix multiplication. 

Let ttij denote the ijth cofactor of M, that is, the determinant of the (n — 1) x (n — 1) matrix 
obtained from M by deleting its ith row and jth column. Let A = Adj(M) denote the adjugate 
matrix of M, that is, the matrix whose ijth entry is {—iy~^^aji. It is a standard consequence of 
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Laplace's determinant expansion [3, §4.111] that MA = (det M)In- Finally, let denote the jth 
column of A. Note that Ma.j = (det M)ej, since both sides are the jth column of (det M)/„. 

Lemma 3. Fix an integer n > 3. Given ann x n matrix M, let aij denote the ijth cofactor of M. 
Also let Z denote the (n — 2) x (n — 2) matrix obtained from M by deleting the first two rows and 
first two columns, so that 

I mil 'IT112 * \ 

(2) 





^ mil 


mi2 


* 




m2i 




* 




\ * 


* 


Z 



Then 011022 ~ 012^^21 — 

(det M) (det Z). 

It is important to note that when det Z ^ 0, the cofactor an is a linear polynomial in the 
variable 777,22 with leading coefficient det Z, while the cofactor 022 is a Unear polynomial in mu 
with leading coefficient det Z (and similarly for the pair 0,12 and 021). For example, when n = 3 
the determinant of the 1x1 matrix Z is simply the lower-right entry 77733 of M\ the identity in 
question is thus 

(7?7ii7?733 - 7?7i37?73i)(7?722"733 - 7?723"732) 

- (7771277733 - 771i377l32) (77l2l77733 - 77723"73i) = 77733 det M. (3) 

For any given dimension n, the assertion of Lemma 3 is simply some polynomial identity that can 
be checked directly; however, a proof that works for all n at once requires a bit of cunning. 

Proof. Define a matrix 

S = (ai a2 63 • • • e„) = 



ail 


-^21 





-ai2 


022 





* 


* 


In-2 



where * represents irrelevant entries. Since B is in lower-triangular block form, its determinant 

det S = det I "^^^ ) • det In-2 = (011O22 - 012O21) 

\~012 022 / 

is easy to evaluate. Moreover, 

MB = (Mai Mai2 Me2, ■ ■ ■ Me^) 

= ((detM)ei (detM)e2 1313 •• • m„) = 



detM 





* 





detM 


* 








z 



Since MB is in upper-triangular block form, its determinant dct(Mi?) = (det M)^(det Z) is also 
easy to evaluate. Using the identity det M • det B = det(M_B), we conclude 

(detM)(oii022 - O12O21) = (detM)^(detZ). 

Both sides of this last identity are polynomial functions of the variables rriij representing the 
entries of M . The factor det M on both sides is a nonzero polynomial, and hence it can be canceled 
to obtain (det M)(det Z) = O11O22 — O12O21 as desired. □ 
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Remark. This proof generalises readily to a similar statement for larger minors of the adjugate 
matrix A. Muir's classic treatise on determinants [10, Ch. VI, §175] includes this generalization of 
Lemma 3 in a chapter wholly devoted to compound determinants (that is, determinants of matrices 
whose elements are themselves determinants). The same result can also be found in Scott's refer- 
ence of equally old vintage [13, p. 62], which has been made freely available online by the Cornell 
University Library Historical Math collection. 

We are now ready to prove Lemma 1. We proceed by induction on n, establishing both n — 2 

and n = 3 as base cases. 

Base case n = 2: The determinant of (™^| ™^^) equals precisely when miim22 = ^12^21. By 
Lemma 2, there are k"^^^ solutions to this equation with the variables mn, mu, ^21, and m22 
all in 2fe. This immediately shows that the number of matrices in Al^l^) is <^e k'^'^^ as claimed. 
Since l/\M.2{k)\ <^ k~^, we see that the probability of a randomly chosen matrix from A42{k) 
being singular is \M'^^{k)\/\M2{k)\ k'^^'. 

Base case n = 3: We first estimate the number of matrices in Ml{k) whose lower right-hand 
entry mss is nonzero. Fix the five entries in the last row and last column of M, with 777,33 7^ 0; there 
are a total of 2k{2k + 1)^ <^ k^ possibiUties. Using the identity (3), we see that if det M — then 
we must have 

(miim33 - mi3m3i)(m22m33 - 77723^32) = (mi2r?733 - 7771377132) (m2im33 - 77l23"^3l)- 

This equation is of the form Li (77711)^2(77122) — Lz{mi2)Li{m2i), where the Li are non-constant 
linear polynomials whose coefficients are at most k"^ in absolute value. (Note that we have used the 
fact that 777,33 7^ in asserting that the Lj are non-constant.) Applying Lemma 2 with a = 2, we 
see that there are <^£; A;^"*"^ solutions to this equation with 777,11, rni2, ^11,21, and 77722 all in J^. This 
shows that there are k'^'^^ matrices in A4'^{k) whose lower right-hand entry 77133 is nonzero. 

If any of the entries in the last row of M is nonzero, then we can permute the columns of M 
to bring that entry into the lower right-hand position; each such resulting matrix corresponds to 
at most three matrices in Ml{k), and so there are still -C^ matrices in Ml{k) that have any 
nonzero entry in the last row. Finally, any matrix whose last row consists of all zeros is certainly 
in Ad^ik), but there are only {2k + 1)^ <^ k^ such matrices. We conclude that the total number of 
matrices in M.'^{k) is A;^+^, so that the probability of a randomly chosen matrix from A43{k) 
being singular is |Al3(/c)|/|A^3(A;)| k~'^~^'^ as claimed. 

Inductive step for 77 > 4: Write a matrix M e M.n{k) in the form (2). Some such matrices will 
have det Z = 0; however, by the induction hypothesis for 77 — 2, the probability that this occurs is 
•Cn,e A;"^^"^ (independent of the entries outside Z), which is an allowably small probability. 

Otherwise, fix values in Xk for the — 4 entries other than m,ii, rrhri, m,2i, and 7^22 such that 
det Z 7^ 0. It suffices to show that conditioning on any such fixed values, the probability that M is 
singular, as 77711, ^12, 77721, and 77122 range over Xk, is <^n,e k~'^'^^. 

By Lemma 3, we see that detM = is equivalent to aiia22 = ai2a2i- Recall that On is a 
linear polynomial in the variable 7^722 with leading coefficient det Z, while the cofactor 022 is a 
linear polynomial in 77711 with leading coefficient det Z (and similarly for the pair 012 and 021). 
Moreover, the coefficients of these linear forms are sums of at most (77,-1)! products of 77 — 1 
entries at a time from M, hence are k"'~^ in size. We may thus apply Lemma 2 with a = 77 — 1 
to see that the probability of 011022 = 012021 is <^n,e k~'^~^^, as desired. □ 

Having established a suitably strong upper bound for (/c) |, we can also bound the cardinality 
of M.n{k) for a fixed A e Z, using the fact that M — XIn will be a singular matrix. Notice that A 
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can be as large as nk if we take M to be the n x n matrix with all entries equal to k, or as small as 
—nk if we take M to be the n x n matrix with all entries equal to —k. It is not hard to show that 
these are the extreme cases for integer eigenvalues, and in fact even more is true: 



Lemma 4. IfM e M.n{k), then every complex eigenvalue of M is at most nk in modulus. 

Proof. Let A be any eigenvalue of M, and let v be a corresponding eigenvector, scaled so that 
maxi<j<„ |vi| = 1 (this is possible since v ^ 0). Then 



|A| = max |(Av)i| = max |(Mv)i| = max 

l<i<n l<i<n l<j<n 



k=l 



Since each entry of M is at most k in absolute value, and each coordinate of v is at most 1 in 
absolute value, we deduce that 



A < max > k^ifcVfc < max > k = nk. 

fc=l k=l 



□ 



Remark. If we use the notation D{z,r)X.o denote the disk of radius r around the complex number 
z, then Lemma 4 is the statement that every eigenvalue of M must lie in D{0,nk). We remark 
that this statement is a weaker form of Gershgorin's "circle theorem" [2], which says that all of the 
eigenvalues of M must lie in the union of the disks 

-Ofmii, ^ DI 17122, ^ ki2ilj, Dirrinn, ^ l"^nj|j- 

l<j<n ^ ^ l<j<n ^ ^ l<J<n ^ 

In fact the proof of Gershgorin's theorem is very similar to the proof of Lemma 4, except that 
one begins with |A — mjj| on the left-hand side rather than |A|, and the other entries of M are left 
explicit in the final inequality rather than estimated by k as above. 

We now have everything we need to prove the Theorem stated earlier: given any integer n > 2 
and any real number £ > 0, the probability that a randomly chosen matrix in M.n{k) has an integer 
eigenvalue is <^n,e l/k^~^. By Lemma 4, any such integer eigenvalue A is at most nk in absolute 
value. For each individual A, we observe that if M e Mn{k) has eigenvalue A, then M — XIn is 
a singular matrix with integer entries which are bounded in absolute value hy k + | A| < {n + l)k. 
Therefore every matrix in M.^{k) is contained in the set 



{M + A7„: MeM°((n + l)A;)}. 



By Lemma 1, the cardinality of this set is <^n,e {{n + 1)A;)" "^"'"^ <^n,e ^" for any fixed 
A. Summing over all values of A between —nk and nk (admittedly, some matrices are counted 
multiple times, but the upper bound thus obtained is still valid), we conclude that the total number 
of matrices in AA^{k) is ^n,£ ^" ^^^^ . In other words, the probability that a matrix in M.n{k) has 
an integer eigenvalue is \M^{k)\/\Mn{k)\ <^n,e l/k^~^, as desired. 
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3. More exact results for 2x2 matrices 



In the special case n = 2, we can sharpen Lemma 1 and the Theorem considerably. The 2x2 
case is particularly nice: since the trace of a matrix with integer entries is itself an integer, it 
follows that if one eigenvalue is an integer then both are. Consequently there is no distinction 
between being diagonalizable over Q and belonging to M^{k). 

Although establishing these sharper results uses only standard techniques from analytic number 
theory, the computations are lengthy and require some tedious case-by-case considerations. There- 
fore we will content ourselves with simply giving the formulas in this section (proofs will appear in 
a subsequent paper [8]). Using the notation f{k) ~ g{k), which means that limfc^oo f{k)/g{k) = 
1, we can state the two formulas in the following way: The probability that a matrix m M2{k) is 
singular is asymptotically 

\Ml{k)\ 6 ^logA: , 



(5) 



\M2{k)\ 7r2 V k^ 

while the probability that a matrix in M.2{k) has an integer eigenvalue is asymptotically 
\Ml{k)\ /7v^ + 4 + 31og(V2 + l)\ /logA;^ 



\M2{k)\ \ 37r2 \ k 



The orders of magnitude (log k) / fc^ and (log k)/k are sharpenings of the orders of magnitude 
^-2-i-£ Lemma 1 and k~^^^ in the Theorem, respectively. 

A consequence of the asymptotic formulas (4) and (5) is that if we choose, uniformly at random, 
only matrices from M.2{k) with integer eigenvalues, then the probability that such a matrix is 
singular is asymptotic to Aa/k, where we have defined the constant 

a = ^ ^ ^ = 0.272008 .... (6) 

14A^ + 8 + 61og(A^+l) 

In other words, the normalized quantity k\M.2{k)\/\M.^{k) \ converges to Aa as k tends to infinity. 
This convergence turns out to be a special case of a more general phenomenon: if we hold A e Z 
fixed and let k tend to infinity, each normalized quantity k\Ai2{k)\/\M^{k) \ converges to the same 
constant Aa. 

However, an interesting picture emerges if instead we rescale A along with k, by thinking of 
A as the nearest integer io 5k with 5 a fixed real number. In fact, there is a continuous function 
U^{S) such that k\M2{k)\/\M^{k)\ ~ U^{X/k) as k and A tend to infinity proportionally to 
each other. The graph of U^{6) is the solid line in Figure 1 below; the exact function is given by 
U'^{5) = aV{\6\), where a is the constant defined in equation (6) and 

^U-2S-S^ + 52log(l + 5) + 2{S - 1) log \5-l\, if < 5 < ^/2, 
\5^-25-\og{5-l)-{5-lf\og{5-l), ifV2<5<2. 

(Note that there is no need to consider V{5) for values of \5\ greater than 2, since all eigenvalues 
of matrices in A42ik) are at most 2k in modulus by Lemma 4.) 

Intuitively, we can think the graph of U^{5) as follows. For each positive integer k, consider 
the histogram of eigenvalues from A1|(A;), vertically normalized by a factor of l/|A^f (A;)|. The 
total area under this histogram is exactly 2, since each matrix contributes exactly two eigenvalues. 
If we then scale the horizontal axis by a factor of 1/ A; and the vertical axis by a corresponding k, 
the total area remains equal to 2, while the horizontal extent of the histogram lies in the interval 
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Figure 1 . Limiting distributions of real and integer eigenvalues for A42{k) 



[—2, 2] . There is one such rescaled histogram for every positive integer k; as k tends to infinity, the 
rescaled histograms converge pointwise to the limiting curve U^{S). 

(The astute reader will notice that we have ignored the fact that matrices with repeated eigenval- 
ues occur only once in (^) but contribute two eigenvalues to the histogram. This effect turns 
out to be negligible: one can show by an argument similar to Lemma 2 that the number of matrices 
in Ai2ik) with repeated eigenvalue A is a vanishingly small fraction of \M2ik)\)- 

For comparison, we can perform the exact same limiting process with the much larger sub- 
set of M.2{k) of matrices having real eigenvalues, which we naturally denote A4f{k). In [4], 
Hetzel et al showed that the probability that a matrix in M.2{k) has real eigenvalues, namely 
(A;)|/|A^2(^)|, converges to 49/72 as k tends to infinity. They observed that this probability 
can be realized as a Riemann sum for the indicator function of the set 



{(a,6,c,d) e M}-. |a|, |c|, \d\ < 1, {a-df + Ahc > 0}, 



where the last inequality is precisely the condition for the matrix ^) to have real eigenvalues. 

If we likewise plot the histogram of eigenvalues from TVff (A:), normalized to have area 2 as 
before, and scale horizontally by l/Zc and vertically by k, we again get a limiting curve U^{6) 
(which bounds an area of exactly 2 as well). The graph of the function U^{6) is the dashed line 
in Figure 1. We can compute this curve using an integral representation similar to the one used 
to derive the constant 49/72; although this integral is significantly more unwieldy, it eventually 
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yields the exact fomula U^{S) = PW{\S\), where P = 72/49 and 



' (80 + 205 + + 525^ - 1075^) /(144(1 + 5)) 
-{5-7S + SS^){1 - 5) log(l - 5) I VI 
-5(l-52)log(l + 5)/4, 



if < 5 < 1, 



W{S) = < 



5(20 + 105 - 1252 - 35=^)/(16(l + 5)) 
+ (35-l)(5-l)log(5-l)/4 
+ 5(52-i)log(5+l)/4, 



if 1 < 5 < v^, 



5(5-2)(2-65 + 352)/(i6(l + 5)) 
-(5-l)3log(5-l)/4, 



if V2 < 5 < 2. 



It is interesting to note the qualitative differences between U and U . Both are even functions, 
since ^A2{k) is closed under negation, and both functions are differentiable, even at 5 = ±\/2 and 
5 = ±2 (except for the points of infinite slope at 5 = ±1). But the real-eigenvalue distribution 
is bimodal with its maxima at 5 ±0.7503, while the integer-eigenvalue distribution is unimodal 
with its maximum at 5 = 0. So a random 2x2 matrix with integer entries bounded in absolute 
value by A; is more likely to have an eigenvalue near 3fc/4 than an eigenvalue near 0; but if we 
condition on having integer eigenvalues, then becomes more likely than any other eigenvalue. 
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